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We study fluctuating tilt Brownian ratchets based on fractional subdiffusion in sticky viscoelastic 
media characterized by a power law memory kernel. Unlike the normal diffusion case the rectification 
effect vanishes in the adiabatically slow modulation limit and optimizes in a driving frequency range. 
It is shown also that anomalous rectification effect is maximal (stochastic resonance effect) at optimal 
temperature and can exhibit a surprisingly good quality. Moreover, subdiffusive current can flow in 
the counter-intuitive direction upon a change of temperature or driving frequency. The dependence 
of anomalous transport on load exhibits a remarkably simple universality. 

PACS numbers: 05.40.-a, 05.10.Gg, 87.16. Uv 



Such diverse research fields as anomalous diffusion and 
transport Brownian ratchets and Stochastic 

Resonance (SR) [l(| attracted much attention over last 
two decades with a huge research literature produced 
and a number of insightful reviews written which ad- 
dress both fundamental aspects of nonequilibrium statis- 
tical physics and various interdisciplinary applications in 
physics, chemistry, biology, and technology. The existing 
ratchet literature is restricted mostly to normal diffusion 
ratchets. Here, a rectification current can emerge for the 
particles diffusing in some periodic and unbiased on av- 
erage potential due to breaking the symmetry of thermal 
detailed balance by an external time-dependent driving. 
This in turn requires to break some spatiotemporal sym- 
metry [H, for example, the spatial inversion symmetr 
as in Fig. [T]in the case of a fluctuating tilt ratchet [T 
driven by harmonic force, which we consider in this work. 
The emergence of net directed motion in unbiased on av- 
erage systems is a strongly nonequilibrium and nonlinear 
effect, absent e.g. within the linear response approxima- 
tion, or linear Onsager regime of nonequilibrium thermo- 
dynamics. A characteristic feature of any true ratchet is 
its ability to sustain a load, i.e. a force directed against 
the transport direction. The presence of a non-zero stop- 
ping force allows to distinguish the genuine ratchets or 
Brownian motors capable to do useful work from the fu- 
tile ones or pseudo-ratchets 0,0. In any isothermal 
Brownian motion which never ceases the dissipative loss 
of energy is compensated at thermal equilibrium by the 
energy (heat) gain due the thermal noise of environment 
so that on average the classical Brownian particle has ki- 
netic energy ksT/2 per degree of freedom. This ensures 
the absence of a net directed motion and of the total heat 
exchange between the particle and its environment. The 
directed Brownian motion requires an external source of 
energy - a part of it will be put into the directed motion 
and a part dissipated as an excess heat to the environ- 



FIG. 1. Ratchet potential and central Brownian particle cou- 
pled to auxiliary Brownian particles modeling viscoelastic en- 
vironment. 



ment. The thermal noise plays in fact a constructive 
role here. First, as a sort of "lubricant" to smooth the 
friction and also providing thermal energy fluctuations 
allowing to overcome potential barriers met on the parti- 
cle's pathway. Restricting to the classical world, without 
noise the Brownian particle would remains localized in 
a potential well, starting there having subthreshold en- 
ergy and for a weak external driving. Therefore, in such 
a setup one generally expects that the rectification cur- 
rent response to a subthreshold driving will increase with 
the noise intensity. However, for a very strong noise the 
potential barriers cease to matter and one expects that 
rectification effect due to a spatial asymmetry of poten- 
tial will vanish. Therefore, there should exist an opti- 
mal thermal noise intensity, or temperature which typies 
Stochastic Resonance (SR) in broad sense. In the focus 
of this Letter are the ratchet and Stochastic Resonance 
effects in subdiffusive transport, i.e. mean displacement 
grows sublinearly, (Sx(t)) cx t a , with < a < 1, in a 
viscoelastic environment. 

The very existence of subdiffusive ratchet transport is 
not obvious and can be questioned. For example, within 
a continuos time random walk (CTRW) mechanism of 
subdiffusion featured by divergent mean residence times 
(MRTs) 0,i the current res pon se to external periodic 
driving is asymptotically zero [11| . This clearly prohibits 
any asymptotic rectification effect for such fluctuating tilt 
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ratchets. However, a subdiffusive rocked ratchet based on 
the fractional Brownian motion (fBm) [l2[ does exist [l3[ , 
and we shall explain its unusual properties, showning in 
particular that the resonance nonadiabatic character of 
this anomalous ratchet effect is of the SR origin, and it 
is a genuine ratchet effect. 

Theory. The Fractional Brownian motion can emerge 
as solution of the Generalized Langevin Equation (GLE) 
for a Brownian particle with mass m, 

mx+ [ <r)(t - t')x(t')dt' = f(x,t)+£(t), (1) 



where f(x,t) = —dV(x,t)/dx is a deterministic force, 
£(f) is zero- mean and Gaussian-distributed thermal noise 
and 77(f) is the frictional memory kernel related to noise 
by the fluctuation-dissipation relation (FDR) 



(mm) = k B Trj(\t-t'\) 



(2) 



The fBm emerges in the overdamped limit, m — > 0, of a 
force-free motion, / — ¥ 0, for a power-law frictional ker- 
nel r](t) = ri a t~ a /r(l - a), with < a < 1 [T(x) is stan- 
dard gamma- function] , and the FDR-related noise £[f), 
which is termed the fractional Gaussian noise (fGN) [12j . 
The corresponding GLE is also termed the fractional 
Langevin equation (FLE) The GLE can be de- 

rived from a standard Hamiltonian model of the Brown- 
ian motion based on the coupling of Brownian particle to 
a thermal bath of harmonic oscillators modeling the envi- 
ronment 15j. The FLE model corresponds to sub-Ohmic 



thermal bath characterized by the spectral bath density 
J(w) oc r; a u; a [l!|. Such a modeling requires a dense spec- 
trum of the thermal bath, or quasi-infinite number of os- 
cillators. Alternatively, one can model the environment 
by a finite set of co-moving auxiliary Brownian parti- 
cles, cf. Fig. [TJ with masses rm coupled elastically with 
constants ki to the central Brownian particle and expe- 
riencing the viscous Stokes friction (with frictional con- 
stants r/i) and uncorrelated, (£,i(t)£j(t')) — 5ij5{t — f'), 
white- noise thermal Gaussian forces y/2r]ikBT^i(t) [20|. 
Considering the overdamped limit for auxiliary particles, 
m, — > 0, this yields upon introduction of fluctuating vis- 
coelastic forces, Ui = —ki{x — Xi) [20j ]: 



«i = —kiV - ViU t + \/2vikik B T£i(t) , (3) 

where i>i = ki/rii are the relaxation rates of viscoelastic 
forces. The last equation for Ui has the form of relax- 
ation equation for elastic force introduced by Maxwell 
in his macroscopic theory of viscoelasticity 2l| which is 
augmented by the corresponding Langevin force in ac- 
cordance with FDR. Such a description was introduced 
m Refs. 0, S3 to model anomalous Brownian mo- 
tion in complex viscoelastic media within a generalized 



Maxwell model. Indeed, by choosing initial itj(0) as in- 
dependent zero-mean Gaussian variables with variance 
(u|(0)) = kiksT and excluding the dynamics of aux- 
iliary variables Uj, the GLE (p}, ([2]) follows immedi- 
ately with the memory kernel 77(f) — J2i=i kie~ Uit . For 
N = 1, an earlier Markovian embedding of GLE with 
exponentially decaying memory kernel [23f is readily re- 
produced. Furthermore, by choosing V{ — i/o/b 1 , ki = 
Caifyrjavf /r(l — a), where b is a scaling parameter 
and C a (b) is a fitting dimensionless constant, the above 
power-law kernel 77(f) can be well approximated over 
about r — N log 10 6 — 2 temporal decades between two 
time cutoffs, 77 = v^ 1 < t < Th = Tib N ~ l . Similar scal- 
ing and approximation are well known in the anomalous 
relaxation theory 0, [24| • Physically, vq corresponds to 
a high-frequency cutoff in J(w), or the largest medium's 
frequency, and v^jb 1 ^^ 1 to the slowest medium's mode, 
which arc always present in reality. For the case a = 0.5 
studied numerically in this work, the choice b = 10, 
Co. 5 = 1.3, N = 12, and i/q = 100 provides an excel- 
lent approximation to the FLE dynamics over at least 
ten time decades, until f max = 10 s . This is checked 
[l3l . 22 . Eil ] by comparison with the exact solution for 



the position variance obtained within and GLE and 
FLE 14, Il7| in the force- free case. We scale time in 



units of the (anomalous) velocity relaxation constant 
t v = (m/77 a ) 1 /( 2_Q ). It is assumed to be temperature- 
independent in accordance with the Hamiltonian model 
TBI lij . This is a standard assumption done also in other 
toy ratchet models @, Q . 

Stochastic dynamics is studied in a driven ratchet po- 
tential, V(x, f) = U(x) — Axcos(Ctt) + fox, where 

U(x) = -Uo[sin{2nx/L) + 0.25 sin^Tnr/L)] (4) 

is a spatially asymmetric periodic potential with ampli- 
tude Uq and period L, A and £1 are the amplitude and 
frequency of the periodic forcing, and fo is a load. The 
distance is scaled in the units of L, the energy in units 
of E = ml? /t%, force in F = E/L and temperature in 
T = E/ks- The role of the inertial effects can be char- 
acterized by the dimensionless parameter r v — l/(uJbT v ), 
where uj h = (2tt / L)(3 3/8 /2 1 ^)^U /m is the bottom and 
(imaginary) top frequency of the considered potential for 
A = 0, fo = 0. The inertial effects can only be negligible 
for r„ < 1 and not too small a [HI]. The corresponding 
borderline value of Uq corresponding to r v = 1 is in the 
dimensionless units Uq w 0.0157. For the simulations 
done in this work the inertial effects are very essential. 
This might seem paradoxical since in the Markovian ap- 
proximation x(t') —> x(t) to Eq. (pj the effective Marko- 
vian friction r] e s(t) — J Q r/(t')dt' oc ri a t l ~ a increases to 
infinity in the course of time. Such a Markovian approxi- 
mation is, however, not affordable for the considered vis- 
coelastic dynamics. Fig. Q] illustrates this point: A part 
of auxiliary particles is strongly coupled to the Brown- 
ian particle. They are co-moving being strongly corre- 
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• U =0.25, n=0.1,AU=0.09,n=2.0 

♦ U =0.645,n=0.1,AU=0.5, n=2.5 
oU =0.645,n=1.0 

AU =1.0, £}=0.1,AU=0.95,n=2.61 
AU =1.0, n=l-0,AU=0.6, n=:2.38 




FIG. 2. (Color online) Anomalous current (subvelocity v a ) as 
a function of temperature for different values Uo and fi. 



lated, reminding a polaron-like picture for quantum par- 
ticles in polar media. However, weaker coupled, stronger 
damped, and much slower particles cannot follow imme- 
diately. They create an elastic force pulling the central 
particle back and retarding its overall motion. This pro- 
hibits any Markovian approximation on the level of the 
reduced (x, v) dynamics as the corresponding slow hid- 
den dynamics cannot be adiabatically eliminated. And 
nevertheless, a highly dimensional Markovian approxi- 
mation with N extra dimensions for overdamped auxil- 
iary particles works remarkably well. Of course, given a 
finite iV the trully asymptotical dynamics becomes nor- 
mal for t ^> t max = b 1 ^^ 1 /^. However, i max grows ex- 
ponentially fast with N and therefore it can be totally 
irrelevant, as in our simulations, to figure out the cor- 
rect asymptotic FLE dynamics. We take a physical limit 
of very large t to study subdiffusive dynamics with i max 
regarded as infinite on this time scale and define the sub- 
velocity as v a = r(l +a) lim feffl' . Practically, the cor- 

t— >oo 

responding values of v a were calculated by fitting the de- 
pendence (x(t)} with the at a dependence, extracting the 
corresponding a within a last time window of simulations 
done until t s j m = 10 6 <C t max . In all simulations we have 
used n = 10 4 trajectories for ensemble averaging. The 
stochastic Heun algorithm has been implemented with 
double precision on a graphical processor unit (GPU) 
26] which allowed to parallelize and accelerate stochas- 
tic simulations by a factor of about 50 as compare with 
standard CPU computing. The amplitude of the external 
driving was fixed at A — 0.8, which is beyond the low- 
est second order of nonlinear response in A/T that yields 
a rectification subtransport response, while the tempera- 
ture, potential amplitude, and the driving frequency were 
varying. 

Results. First, we fixed the temperature at T — 0.25 
and varied the driving frequency and the potential am- 
plitude, cf. [20l |. Fig. SI. Remarkably, the rectification 



resp onse is absent in the adiabatic limit Q, — > as in 
[1 31 ] . This is in a sharp contrast with the case of normal 
diffusion where the adiabatic current is maximal [H, 0] ■ 
Indeed, the sub-diffusion and subtransport in periodic po- 
tentials turn out to be asymptotically insensitive to the 
potential amplitude Uq [221 ]. This surprising feature is 
due to the influence of slugish dynamics of the medium's 
degrees of freedom which cannot immediately follow to a 
faster moving Brownian particle. They cause ultraslow 
dynamics on the time scale which largely exceeds the 
mean time spent in a potential well. In this respect, the 
medium's dynamics is not influenced directly by an exter- 
nal force. Therefore, for an adiabatically slow driving a 
periodic potential does not play any role in the long time 
limit and the ratchet transport is absent. For a very fast 
driving, the transport is also obviously increasingly sup- 
pressed upon increasing the driving frequency. Therefore, 
there should exist an optimal value of driving frequency 
when the corresponsing subvelocity attains a maximum, 
v a ,max — maxo v n (fl), which is a Stochastic Resonance 
related effect [20j. This maximal value is optimized also 
with the potential amplitude. For example, v a . max is 
larger for Uq = 0.645, than for Uq — 0.5 and Uq = 1.0, 
see in 1201 ] . Fig. SI and the corresponding discussion 
therein. Noteworthy is also the inversion of the current 
direction for sufficiently large Uq and f2. The subtrans- 
port then occurs in the conterintuitive direction, contrary 
to the direction predicted by the slow adiabatic tilt argu- 
mentation. A similar inversion for high-frequency driving 
occurs also in the normal diffusion case 0, H| ■ 

SR-related phenomenon occurs also in the depen- 
dence of subcurrcnt on temperature, cf. Fig. [2] 
An adiabatic driving argumentation predicts v a (T) oc 
T-"exp(-AC//T) with n = 2 and AU = 2U in the 
lowest order of perturbation theory over A/T (quadratic 
response), for a weak driving, AL <C 2ttUq. This is a typ- 
ical SR dependence yielding the maximum versus tem- 
perature at T max = AU/n. We are dealing but with a 
strong nonadiabatic driving and hence use AU and n 
in Fig. [2] as some fitting parameters. The SR origin 
of the maxima for sufficiently small f2 is substantiated 
by matching the forward tilt half-periods with the nu- 
merical mean times (t\j) of the jumps into the forward 
directions, i.e. (ri,/) = w/Q, [l(J, with temperature near 
to the transport optimization as indicated by arrows in 
Fig. [21 The negative current for a larger ft is also op- 
timized with temperature, which does not have but any 
relation to SR or synchronization. Morever, an inversion 
of the subcurrent direction with temperature is detected 
in Fig. [2] (see also the insert therein) for Uq = 0.645 and 
a high-frequency driving fl = 1.0. Similar temperature 
inversions occur also in the normal diffusion case Q. 

Furthermore, the viscoelastic sub-diffusion in periodic 
potentials is asvmptotically also not affected by the po- 
tential height 
slightly affect it 



22j even though a periodic driving can 
Asymptotically, in the considered units 
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FIG. 3. (Color online) Anomalous current (subvelocity v a ) as 
a function of positive load fo at different Uo, T and a small 
driving frequency £1 — 0.1. Insertion shows dependencies v a 
versus negative load fo at T = 0.3 and some high Q — 1.0 at 
different Uo. 



(6x 2 (t)) ~ 2D a t a /T(l + a) with subdivision coefficient 
D a = T. A periodic driving does affect D„. However, 
it is not changed strongly [13|. Therefore, the general- 
ized Peclet number (l~3jj Pe Q := v a L/D a , which measures 
the coherence and quality of transport [27| . can be ap- 
preciable. With lowering the barrier height in Fig. [2] 
the maximum of subvelocity is shifted towards smaller 
temperatures and the corresponding Pe a can exceed the 
value of one signifying thereby that a high-quality subd- 
iffusivc ratchet transport is possible. 

Finally, the dependence of the directed subtransport 
on the load fo in the opposite direction is shown in 
Fig. [3] The existence of the stopping force shows 
clearly that we are dealing with a genuine ratchet ef- 
fect. Given the asymptotic independence of the viscoelas- 
tic subtransport on the potential amplitude Uq in wash- 
board potentials (22j one expects a very simple depen- 
dence v a (fo) = v a (0) — fo (in dimesionless units) to hold 
for small driving frequencies il. Indeed, the numerical 
results are consistent (taking numerical errors into ac- 
count) with this prediction. However, for larger f2 and 
the inverted transport the deviations become more pro- 
nounced. 

In conclusion, with work we put forward a toy model 
for viscoelastic subdiffusive ratchet transport and showed 
that such a transport can be optimized by ambient ther- 
mal noise and frequency of the external driving due to 
Stochastic Resonance effects and can exhibit a surpris- 
ingly good quality. 
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Figure SI: Anomalous current (subvelocity v a ) as a function CI 

of the driving frequency for different Uq at T = 0.25 for zero Figure S2: Mean times of transitions, (tij) and (r 2 ,/), as 
load fo — 0. functions of the driving frequency f2. The values f2 max and 

t) a (!! max ) are indicated by crosses. 



Supplementary Material 

1. Stochastic equations of motion for a central particle 
surrounded by a cloud of elastically bounded Brownian 
particles read 

JV 

mi = f(x, t) - ^2 ki(x - Xi) , 
i=i 

TUiXi = ki(x ~ Xi) - rjiXi + \j2r)ik B T^i(t) . (5) 
This equations are similar to the starting equations of 
motion in the purely dynamical model leading to the 
GLE description. The last one can be obtained by setting 
rji — > 0, i.e. considering a purely dynamical evolution with 
a tremendously large number of oscillators N — ¥ oo with 
a dense spectrum u)i — y/ki/rru- The only nondynami- 
cal assymption in the dynamical model is that initially 
the thermal bath oscillators are canonically distributed 
at the temperature T - like in a standard molecular dy- 
namics setup. We replace them by a finite set of auxiliary 
Brownian particles. 

2. Fig. SI displays the dependence of subvelocity on ex- 
ternal driving frequency Q. for different potential ampli- 
tudes Uq and the fixed driving amplitude and temper- 
ature. In order to clarify the physical origin of maximal 
v a ,m.ax we have calculated the frequency-dependent mean 
times of the transitions to the neighboring potential well, 
(tij), and to the second next potential well, (ra,/). If the 
maximum of v a is due to synchronization, then the trans- 
port should become optimal when the particles advance 
in the trasport direction over one or two potential pe- 
riods during the driving half-period 71/2 = 7r/f2. This 
means that the condition (tij) < 71/2 < ( r 2,/) should 
be obeyed. Indeed, the numerical results displayed in Fig. 
S2 show that this is the case for not too high Uo- How- 
ever, since (tij) and (t%j) are only weakly frequency- 
dependent and thus no frequency entraintment occurs, 
we are dealing with the phenomenon of Stochastic Reso- 
nance rather than stochastic synchronization. Here, the 
response maximizes when an external driving frequency 
fits into an intrinsic frequency of stochastic dynamics. 



